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The CIE 1997 Interim Colour Appearance Model 
(Simple Version), CIECAM97s 

SUMMARY 

TCI -34. Testing Colour Appearance Models, was established to test various 
models for the prediction of the colour appearance of object colours. Later, 
TCl-34's terms of reference were extended to include the recommendation of 
a single colour appearance model to be used until a demonstrably better 
model could be formulated. The committee decided to formulate a single 
model with a simple version for many practical applications and 
comprehensive version for a wide range of viewing conditions and 
phenomena. TCI -34 will continue its work testing this and other models. 
This report summarizes the formulation of the simple version of the CIE 
Interim Colour Appearance Model, CIECAM97s. The extension of this model 
to a comprehensive version, CIECAM97c, will be formulated and published 
in the future. 
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1. INTRODUCTION 

In March 1996 the CIE held an expert symposium on Colour Standards for 
Image Technology in Vienna (CIE, 1996). While the symposium covered 
many image technology concepts for which the CIE might provide guidance 
or standards to assist industry, one of the critical issues was the establishment 
of a general-use colour appearance model. Symposium participants from 
industry recognized the need for a colour appearance model, but requested 
CIE guidance in establishing a single model that could be used throughout the 
industry to promote uniformity of practice and compatibility between various 
components in modern imaging systems. 

1.1 Background 

Colour appearance models are used to extend traditional colorimetry (e.g., CIE 
XYZ and CIELAB) to the prediction of the observed appearance of coloured 
stimuli under a wide variety of viewing conditions. This is accomplished by 
taking into account the tristimulus values of the stimulus, its background, its 
surround, the adapting stimulus, the luminance level, and other factors such 
as cognitive discounting the illuminant. The output of colour appearance 
models includes mathematical correlates for perceptual attributes such as 
brightness, lightness, colourfulness, chroma, saturation, and hue. Colour 
appearance models are used in applications such as digital colour 
reproduction and assessment of the colour rendering properties of light 
sources. 

The requirements for a single model were highlighted and summarized in a 
presentation by Hunt (1996) made at the Vienna symposium. In that 
presentation, Hunt reviewed the status and historical development of 
various models and presented 12 principles for consideration in establishing a 
model. These principles, reproduced below, served as the guiding rules in the 
formulation of CIECAM97s. 
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1. The model should be as comprehensive as possible, so that it can be 
used in a variety of applications; but at this stage, only static states of 
adaptation should be included, because of the great complexity of 
dynamic effects. 

2. The model should cover a wide range of stimulus intensities, from 
very dark object colours to very bright self-luminous colour. This 
means that the dynamic response function must have a maximum, 
and cannot be a simple logarithmic or power function. 

3. The model should cover a wide range of adapting intensities, from 
very low scotopic levels, such as occur in starlight, to very high 
photopic levels, such as occur in sunlight. This means that rod vision 
should be included in the model; but because many applications will be 
such that rod vision is negligible, the model should be usable in a 
mode that does not include rod vision. 

4. The model should cover a wide range of viewing conditions, 
including backgrounds of different luminance factors, and dark, dim, 
and average surrounds. It is necessary to cover the different surrounds 
because of their widespread use in projected and self-luminous 
displays. 

5. For ease of use, the spectral sensitivities of the cones should be a 
linear transformation of the CIE x, y, z or x l0 , y l0 , z, 0 functions, and 
the V f (X) function should be used for the spectral sensitivity of the 
rods. Because scotopic photometric data is often unknown, methods of 
providing approximate scotopic values should be provided. 
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6. The model should be able to provide for any degree of adaptation 
between complete and none, for cognitive factors, and for the Helson- 
Judd effect, as options. 

7. The model should give predictions of hue (both as hue-angle, and as 
hue-quadrature), brightness, lightness, saturation, chroma, and 
colourfulness. 



8. The model should be capable of being operated in a reverse mode. 

9. The model should be no more complicated than is necessary to meet 
the above requirements. 

10. Any simplified version of the model, intended for particular 
applications, should give the same predictions as the complete model 
for some specified set of conditions. 

11. The model should give predictions of colour appearance that are 
not appreciably worse than those given by the model that is best In each 
application. 

12. A version of the model should be available for application to 
unrelated colours (those seen in dark surrounds in isolation from 
other colours). 

The conclusion drawn at the symposium was that the CIE should 
immediately begin work on the formulation of such a model with the goal 
that it be completed prior to the AIC (International Colour Association) 
quadrennial meeting to be held in Kyoto in May, 1997. The CIE decided that 
this task would become part of the TCI -34 terms of reference. This report 
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details the simple version of the model agreed upon by TCI -34 at its meeting 
during the May, 1997 CIE Division 1 meeting in Kyoto. 

There are a number of terms that require precise usage in the colour 
appearance field. All terms in this report follow the definitions of the 
International Lighting Vocabulary (CIE, 1987). Their typical usage in the 
colour appearance field is also further detailed in Fairchild (1998). 

2. THE CIECAM97s MODEL 

It is important to note that the formulation of CIECAM97s builds upon the 
work of many researchers in the colour appearance field. This was an 
important issue in TCl-34's establishment of this model as the best of what is 
currently available. Various aspects of the model can be traced to work of (in 
alphabetical order) Bartleson. Breneman, Fairchild, Estevez, Hunt, Lam. Luo, 
Nayatani, Rigg. Seim. and Valberg among others. Examples of these 
contributions include the Bradford chromatic-adaptation transform (Lam, 
1985; Luo, 1997), a different exponent on the short wavelength response 
(Nayatani et ah, 1982), partial adaptation factors (Fairchild. 1996; Nayatani, 
1997), cone responsivities (Estevez; see Hunt and Pointer. 1985). a hyperbolic 
response function (Seim and Valberg, 1986). redness-greenness and 
yellowness-blueness scales (Hunt. 1994; Nayatani. 1995). surround effects 
(Bartleson and Breneman, 1967). no negative lightness predictions (Nayatani, 
1995, Fairchild. 1996). and a chroma scale (Hunt, 1994). Clearly, CIECAM97s is 
the result of the amalgamation of a wide range of colour appearance research. 

Preliminary testing of CIECAM97s indicates that it performs as well as. if not 
better, than any previously published model for a wide range of experimental 
data (Hunt, 1997). These results will be expanded upon and published in the 
final report of TCI -34. The comprehensive model. CIECAM97c, will be 
derived from the simple model, CIECAM97s. by adding features that include 
(but are not necessarily limited to) the rod response, prediction of the Helson- 
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Judd and Helmholtz-Kohlrausch effects, and options for application to 
unrelated colours. It is also important to note that these models are 
recognized to be empirical models capable of predicting the available visual 
data. It is expected that anticipated scientific insights will, at some time in the 
future, allow the derivation of more theoretically correct models that are also 
capable of predicting experimental results. 

2.1 Input Data 

The model input data are the adapting field luminance in cd/m 2 (normally 
taken to be 20% of the luminance of white in the adapting field). L A , the 
relative tristimulus values of the sample in the source conditions, XYZ, the 
relative tristimulus values of the source white in the source conditions, 
XwY w Z w , and the relative luminance of the source background in the source 
conditions, Y b . Additionally, the constants c, for the impact of surround, N c , a 
chromatic induction factor, F LL , a lightness contrast factor, and F, a factor for 
degree of adaptation, must be selected according to the guidelines in Table 2.1. 
All CIE tristimulus values are obtained using the CIE 1931 Standard 
Colorimetric Observer (2 ) (CIE, 1986). Background is defined as the area 
immediately adjacent to the stimulus of interest and surround is defined as 
the remainder of the visual field. Surround relative luminances of greater 
than or approximately equal to 20% of the scene white are considered average, 
less than 20% are considered dim, and approximately 0% are considered dark. 



Table 2.1. Selection guidelines for parameters used in CIECAM97s. 



Viewing Condition 


c 


N c 


Fa. 


F 


Average Surround, Samples Subtending> 4 ° 


0.69 


1.0 


0.0 


1.0 


Average Surround 


0.69 


1.0 


1.0 


1.0 


Dim Surround 


0.59 


1.1 


1.0 


0.9 


Dark Surround 


0.525 


0.8 


1.0 


0.9 


Cut-Sheet Transparencies (on a viewing box) 


0.41 


0.8 


1.0 


0.9 
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2.2. Chromatic Adaptation 

An initial chromatic adaptation transform is used to go from the source 
viewing conditions to corresponding colours under the equal-energy- 
illuminant reference viewing conditions. First, tristimulus values for both 
the sample and white are normalized and transformed to spectrally- 
sharpened cone responses, illustrated in Fig. 2.1 ., using the transformation 



R(H) 
G(H) 
B(H) 



— R(B) 

G(B) 

B(B) 
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.Tnes) ,0n ' ( ShCd lin6S) Hu ^ Pointer-Estevez transformation 4 (soHd 
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given in Eqs. 2.1 and 2.2. Note that the forward matrix transformation given 
in Eq. 2.2 was applied to the spectral tristimulus values of the CIE 1931 
Standard Colorimetric Observer in order to generate the curves presented in 

Fig. 2.1. 



tri rx/Yi 

IgI = M b Iy/yI (2.1) 

W Lz/yJ 

r 0.8951 0.2664 -0.16141 I" 0.9870 -0.1471 0.16001 

M B = 1-0.7502 1.7135 0.0367 1 M B -1 = I 0.4323 0.5184 0.0493 1(2.2) 

[ 0.0389 -0.0685 1.0296 J (_-00085 0.0400 0.9685 J 



The chromatic-adaptation transform is a modified von Kries-type 
transformation with an exponential nonlinearity on the short-wavelength 
sensitive channel as given in Eqs. 2.3 through 2.6. In addition, the variable D 
is used to specify the degree of adaptation. D is set to 1.0 for complete 
adaptation or discounting the illuminant (as is typically the case for reflecting 
materials). D is set to 0.0 for no adaptation. D takes on intermediate values 
for various degrees of incomplete chromatic adaptation. Equation 2.7 allows 
calculation of such intermediate D values for various luminance levels and 
surround conditions. 



R c =[D(1.0/R W ) + 1-D]R (2.3) 

G C =[D(1.0/G W ) + 1-D]G (2.4) 

B c =[D(l.0/B^) + l-D]B| P (2.5) 

P = (B w /1.0r* (2.6) 

D = F- F/[l + 2(L^ ) +(L 2 A )/300] (2.7) 



If B happens to be negative, then B c is also set to be negative. Similar 
transformations are also made for the source white since they are required in 
later calculations. Various factors must be calculated prior to further 
calculations as shown in Eqs. 2.8 through 2.12. These include a background 
induction factor, n, the background and chromatic brightness induction 
factors, N bb and N cb , and the base exponential nonlinearity, z. 

k = l/(5L A+ l) (2.8) 

F L = 0.2k<(5L A )+0.l(l-k<) 2 (5L A ) 1/3 (2.9) 

n = Y b /Y w (2.10) 

N bb = N eb =0.725(l/n) a2 (2 .n) 

Kin l/Z (2.12) 



z=l + F lin ,/2 



The post-adaptation signals for both the sample and the source white are then 
transformed from the sharpened cone responses to the Hunt-Pointer-Estevez 
cone responses as shown in Eqs. 2.13 and 2.14, and illustrated in Fig. 1, prior to 
application of a nonlinear response compression. 

TR'l Tr c y1 

JG'I-MhMS'Ig.yI (2.13) 

Lb 1 J [b c yJ 

r 0.38971 0.68898 -0.078681 Tl.9102 -1.1121 0.20191 

M H = 1-0.22981 1.18340 0.04641 I = 1 0.3710 0.6291 0.00 1(2.14) 

L 0.00 0.00 1.00 J [ 0.00 0.00 1.00 J 

The post-adaptation cone responses (for both the sample and the white) are 
then calculated using Eqs. 2.15 through 2.17. 



\0.73 



pi 40(F L RV100)° 

a "[(F L RV100r+2] +1 (215) 
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40(F L G'/100) 073 



+ 1 



(2.16) 



, 40(F L B'/100)°" 
a ~[(F L BVlOO) 0 73 + 2] 



+ 1 



(2.17) 



2.3. Appearance Correlates 

Preliminary red-green and yellow-blue opponent dimensions are calculated 
using Eqs. 2.18 and 2.19. 



The CIECAM97s hue angle, h, Is then calculated from a and b using Eq. 2.20. 



Hue quadrature, H, and eccentricity factors, e, are calculated from the 
following unique hue data via linear interpolation between the following 
values for the unique hues: 

Red: h = 20.14, e = 0.8, H « 0 or 400, 

Yellow: h = 90.00, e = 0.7, H = 100, 

Green: h = 164.25, e = 1.0, H = 200, 

Blue: h = 237.53. e = 1.2. H = 300 

Equations 2.21 and 2.22 illustrate calculation of e and H for arbitrary hue 
angles where the quantities subscripted 1 and 2 refer to the unique hues with 
hue angles just below and just above the hue angle of interest. 



a = R' 3 -l2G' ii /ll+B',/11 
b=(l/9XR' 3 +G , a -2B'J 



(2.19) 



(2.18) 



h = tan"'(b/a) 



(2.20) 



e = e 1 +(e 2 -e,)(h-h,)/(h 2 -h 1 ) 



(2-21) 
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„ u 100(h-h.)/e, 

H = H . + *- — i (? 22) 

1 (h-h,)/e, + (h 2 -h)/e 2 K } 



The achromatic response is calculated as shown in Eq. 2.23 for both the 
sample and the white. 

A = [2R 1 . +G' , + (1/20)B 1 . -2.05]N bb (2.23) 

CIECAM97s Lightness, J. is calculated from the achromatic signals of the 
sample, A. and white, A w , using Eq. 2.24. 

J = 100(A/A w f (2.24) 

CIECAM97s Brightness, Q. is calculated from CIECAM97s lightness and the 
achromatic response for the white using Eq. 2.25. 

Q=(1.24/cXJ/100r 7 (A w+ 3)° 9 (2.25) 

Finally, CIECAM97s saturation, s; CIECAM97s chroma, C; and CIECAM97s 
colourfulness, M; are calculated using Eqs. 2.26 through 2.28, respectively. 



50(a' +b 2 ) lOOedO/lSjN.N^ 

R' a +G , a + (21/20)B , a (2 - 26) 

C= 2.44s 069 (J/100) 067n (l.64 -0.29") (2.27) 

M=CF L 0,S (2.28) 



i 1 

3. CONCLUSIONS 

CIE TCI -34 recommends that the CIECAM97s model be evaluated as an 
interim solution to the problem of colour appearance specification. This 
model should help to address industrial needs by providing a single, CIE- 
recognized colour appearance model that represents the committee's best 
effort to bring together the best features of existing models. TCI -34 plans to 
subject this and other models to further testing. It is reasonable to expect that, 
at some future date, a more accurate and/or theoretically-based model might 
be developed. 

3.1. CIECAM97s (Simple Version) Recommended Use 

The CIECAM97s model should be adequate for most practical applications 
requiring use of colour appearance metrics or transformations more 
sophisticated than those provided by the CIELAB colour space. It does not 
allow for predictions of the influence of rod photoreceptors on colour 
appearance, the Helson-Judd effect, the Helmholtz-Kohlrausch effect, or the 
appearance of unrelated colours. A more comprehensive model, such as the 
planned CIECAM97c model should be considered when such phenomena are 
important. 

CIECAM97s provides mathematical scales to correlate with various 
perceptual appearance attributes. As such, it does not explicitly construct a 
colour space. The CIECAM97s lightness, chroma, and hue correlates G,C,h) 
can be used to construct a colour space by considering them as cylindrical 
coordinates as is done in the CIELAB colour space with L*. C* ab , and h ab . 
Alternatively, a brightness-colourfulness space could be constructed using 
CIECAM97s Q, M, and h as cylindrical coordinates. If a rectangular space is 
required, one can be constructed using the normal means for cylindrical-to- 
rectangular coordinate transformations (i.e., J, Ccos(h), and Csin(h) or Q, 
Mcos(h), and Msin(h) could be used as rectangular coordinates.) 
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The question of colour-difference specification is often closely linked to that 
of colour appearance modeling. At this time, the CIECAM97s model has not 
been evaluated for use as a colour-difference space. While it would be 
possible to calculate a Euclidean colour-difference metric in either of the 
rectangular colour spaces described above, such a practice has not yet been 
shown to be better or worse than the current recommendation for CIE94 
colour differences (CIE, 1995). 



13 

4. APPENDIX 

4.1. Numerical Examples 

Example calculations using the CIECAM97s are given for four samples in 
Table 4.1. Note that the magnitudes of the lightness and chroma scales in 
CIECAM97s are comparable to those of CIELAB. All calculations were 
performed using double precision. Slightly different results might be 
obtained with single precision calculations. (The hue angle calculation for 
near neutral colours (e.g., Case 1) is a particular example in which large 
numerical differences can arise that have little visual meaning since the hue 
of an achromatic sample is truly undefined.) A Microsoft Excel spreadsheet 
with these example calculations for CIECAM97s can be found at 
<http://www.cis.rit.edu/people/faculty/fairchild/CAM.html>. 

Table 4.1. Example calculations using CIECAM97s for four samples. 





Case J 


Case 2 


Case 3 


Case 4 


X 


19.01 


57.06 


3.53 


19.01 


Y 


20.00 


43.06 


6.56 


20.00 


Z 


21.78 


31.96 


2.14 


21.78 




95.05 


95.05 


109.85 


109.85 




100.00 


100.00 


100.00 


100.00 


z w 


108.88 


108.88 


35.58 


35.58 


L A (cd/m 2 ) 


318.31 


31.83 


318.31 


31.83 


F 


1.0 


1.0 


1.0 


1.0 


D 


0.997 


0.890 


0.997 


0.890 


Y b 


20.0 


20.0 


20.0 


20.0 


c 


0.69 


0.69 


0.69 


0.69 


N c 


1.0 


1.0 


1.0 


1.0 




1.0 


1.0 


1.0 


1.0 


k 


0.0006 


0.0062 


0.0006 


0.0062 




1.17 


0.54 


1.17 


0.54 
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U.cU 




0.20 


0.20 


N 

1N bb 


l.UU 


1.00 


l.OO 


1.00 


M 

1N cb 


i nr. 


l.OO 


l.OO 


1.00 


Z, 




1.45 


1.45 


1.45 


i\ 


u.y4 


1.33 


0.70 


0.94 


P 


1.U4 


0.75 


1.32 


1.04 


D 


1.09 


0.75 


0.29 


1.09 




0.94 


0.94 


1. 19 


1.19 




1.04 


1.04 


0.90 


0.90 


B w 


1.09 


1. 09 


0.34 


0.34 


p 


1.01 


l.Ol 


0.91 


0.91 


D 


1.00 


l.4l 


0.58 


0.81 


G e 


1.00 


0.72 


1.46 


1.14 


13 

B c 


1.00 


0.69 


0.86 


2.70 


K CW 


1.00 


0.99 


1.00 


1.02 


G cw 


1.00 


l.OO 


1.00 


0.99 




1.00 


l.Ol 


1.00 


0.93 


K 


20.0 


51.2 


5.58 


18.9 


G 


20.0 


39.3 


7.70 


23.0 


B 


20.0 


29.5 


5.80 


53.0 


K w 


100.0 


99.7 


100.0 


101.0 


G w 


100.0 


100.2 


100.0 


99.4 


B w 


100.0 


lOl.O 


99.8 


93.3 


R a 


6.90 


7.56 


3.55 


4.46 


° a 


D.90 


6.57 


4.17 


4.94 


B a 


6.90 


5.64 


3.62 


7.70 




15.4 


10.7 


15.4 


10.7 


G' 

aw 


15.4 


10.7 


15.4 


10.7 


B' 

aw 


15.4 


10.7 


15.3 


10.3 


a 


-0.0005 


0.90 


-0.67 


-0.23 
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1- 

u 


-U.UUU4 


0.32 


0.05 


-0.67 


n 


^1 y.4 


19.35 


175.4 


250.8 


T T 
H 




399 


218 


307 


H c (Red) 


u 


99 


0 


7 


T T / v n > 

ri c (Yellow; 


U 


0 


0 


0 


ri c (Green) 


30 


0 


82 


0 


H c (Blue) 


70 


1 


18 


93 


e 


1.15 


0.80 


1.03 


1.16 


A 


18.99 


19.92 


9.40 


12.19 


A 

w 


44.80 


30 54 


44 80 




J 


42.44 


65.27 


21.04 


39.88 


Q 


32.86 


31.88 


20.53 


22.96 


s 


0.14 


146.98 


232.16 


180.56 


C 


0.47 


61.97 


72.99 


66.85 


M 


0.49 


56.52 


74.70 


60.98 
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4.2. Inverting the CIECAM97s Model 

Steps for using the CIECAM97s model In the reverse direction for 
corresponding-colours calculations or colour-reproduction applications 
follow. 

Starting Data: 

QorJ, MorC. Horh 

A w , n. z, F L> N bb , N cb Obtained Using Forward Model 
Surround Parameters: F, c, F LL , N c 
Luminance Level Parameters: L A , D 
Unique Hue Data: 

Red: h = 20.14. e = 0.8 

Yellow: h = 90.00, e = 0.7 

Green: h = 164.25, e= 1.0 

Blue: h = 237.53, e = 1.2 

(1) From Q Obtain J (if necessary) 

J = 100(Qc/ 1.24) 1/067 /(A w + 3)° 9/0 67 ( 4 . 10 

(2) From J Obtain A 

A=(J/100)"=A W (4.2) 

(3) Using H. Determine h, , h 2 , e, . e 2 (if h is not available) 

e, and h, are the values of e and h for the unique hue having the 
nearest lower value of h and e 2 and h 2 are the values of e and h for the 
unique hue having the nearest higher value of h. 

(4) Calculate h (if necessary) 

h = [(H-H ) )(h ) /e I -h 2 /e 2 )-100h,/e 1 )4(H-H 1 )(l/e I -l/e 2 )-100/e 1 ] 

(4.3) 
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H, is 0, 100. 200, or 300 according to whether red, yellow, green, 
or blue is the hue having the nearest lower value of h. 

(5) Calculate e 

e = e, + (e 2 -e 1 )(h-h,)/(h 2 -h,) (4.4) 
e ] and h, are the values of e and h for the unique hue having the 
nearest lower value of h and e 2 and h 2 are the values of e and h 
for the unique hue having the nearest higher value of h. 

(6) Calculate C (if necessary) 

C=M/F L 015 (4.5) 

(7) Calculate s 

s=c ./o.69 /[2.44Q/100) 067n (1.64-0.29 n )] ,/069 (4.6) 

(8) Calculate a and b 
a=s(A/N^ + 2.05)/{[l+(tanty 

(47) 

In calculating [l +(tan h) 2 ]" 2 the result is taken as: 
positive for 0 °< h < 90 ° 
negative for 90 °< h < 270 0 
positive for 270 °< h < 360 °. 



b = a(tanh) (4.8) 

(9) Calculate R' a , G' a , and 

R' a = (20/61)(A/N bb + 2.05) + (41/61)(ll/23)a+ (288/61) (1 /23)b (4.9) 
G' a = (20/61)(A/N bb +2.05)-(81/61)(ll/23)a-(261/61)(l/23)b (4.10) 
B' a = (20/61)(A/N bb + 2.05)-(20/61)(ll/23)a-(20/61)(315/23)b (4.11) 
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(10) Calculate R', G'. and B' 

R = 100[(2R' 4 -2)/(41-R' a )] 1/ 
G = 100[(2G' a -2)/(41-G , a )] l/073 
B , = 100[(2B , 1 -2)/(41-B' a )] 1/0n 

If R' a -1 < 0 use: 

R' = -100[(2- 2R' a )/(39 + R' , )]' /073 
and similarly for the G' and B' equations. 



(4.12) 
(4.13) 
(4.14) 

(4.15) 



(11) Calculate R C Y, G C Y, and B C Y 



R C Y 




RVF L 


G C Y 


= M B M" H ' 


GVF L 


B C Y 




BVF L 



(4.16) 



(12) Calculate Y c 



Y c =0.43231R C Y + 0.51836G C Y + 0.04929B C Y 



(4.17) 



(13) Calculate (Y/Y C )R. (Y/YJG, and (Y/YJ'^B 

(Y/Y c )R = (Y/Y e )R c /[D(l/RJ + l-D] 
(Y/Y C )G=(Y/Y C )G C /[D(1/GJ+ 1-D] 
(Y/ Y c ) ,/p B = [|(Y/Y C )B C |]' /P /[D(l /B * ) + 1 - D] 
If (Y/Y e )B c <0.0 then (Y/YJ'^B is also set to be negative. 



i/ P 



(4.18) 
(4.19) 
(4.20) 



(14) Calculate Y' 



Y' = 0.43231YR+0.51836YG + 0.04929(Y/Y C )" P BY 



(4.21) 
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(15) Calculate X". Y" and Z" 



X' ' 




Y C (Y/Y C )R 




Y" 


=m;' 


Y C (Y/Y C )G 


(4.22) 


Z' ' 




Y c (Y/Y c )' /p B/(Y'/Y c ) (,/p -" 





Note: X'\ Y", and Z" are equal to the desired X, Y, and Z to a very close 
approximation. This is because Y 1 differs from Y since (Y/Y c ) 1/p BY c is used 
instead of YB. However this is multiplied by 0.04929 so the difference is 
small. 
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